Rotational Fluxons of Bose-Einstein Condensates in Coplanar Double-Ring Traps 
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Rotational analogs to magnetic fluxons in conventional Josephson junctions are predicted to emerge in 
the ground state of rotating tunnel-coupled annular Bose-Einstein condensates (BECs). Such topological 
condensate-phase structures can be manipulated by external potentials. We determine conditions for observ- 
ing macroscopic quantum tunneling of a fluxon. Rotational fluxons in double-ring BECs can be created, ma- 
nipulated, and controlled by external potential in different ways than possible in the solid state system, thus 
rendering them a promising new candidate system for studying and utilizing quantum properties of collective 
many-particle degrees of freedom. 
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Remarkable experimental advances have made it possible 
to engineer cold atom systems to represent landmark mod- 
els from completely different fields of physics. Examples 
include quantum phase transitions Jl]] and the Josephson ef- 
fect H. Besides intriguing nonlinear dynamics, the Joseph- 
son effect shows macroscopic quantum phenomena with ex- 
citing prospects for applications J^. Long Josephson junc- 
tions were used, e.g., to trap and study magnetic flux quanta, 
and the macroscopic quantum tunneling of such fluxons was 
observed J3]. Here we report the existence of topologi- 
cal condensate-phase structure equivalent to superconducting 
fluxons in rotating BECs that are confined in two concen- 
tric ring-shaped traps. The BECs in the individual rings are 
coupled by tunneling through a potential barrier. The rota- 
tional fluxons can be understood as vortices that have entered 
the tunnel barrier. They show intriguing dynamical behav- 
ior and macroscopic quantum properties. Easier accessibility 
and more straightforward means of manipulation than possi- 
ble in conventional Josephson junctions make rotational flux- 
ons in tunnel-coupled BECs attractive for investigating funda- 
mental problems ranging from models for cosmological evo- 
lution 15[] topossibilities for realizing quantum information 
processing jfj. 

Recent successful efforts to create annular trapping geome- 
tries 171 l8D and the routine use of trap rotation to simulate the 
effect a magnetic field has on charged particles J9(] have mo- 
tivated our present theoretical work. Unlike in the previously 
considered cases of vertically separated double-ring traps II 1 Oil 
or coupled elongated BECs jl ll Il2ll . we predict rotational 
fluxons to occur in the ground states of the proposed system. 
In contrast to unconfined vortices in harmonically trapped 
BECs, rotational fluxons are confined to the tunnel barrier 
region between the coupled rings for energetic reasons and 
thus take on properties of topological solitons [13]. Preparing 
ground state solitons by cooling opens unprecedented oppor- 
tunities for precision experiments on classical and quantum 
soliton dynamics. The phase structures are analogous to mag- 
netic flux quanta occurring in a superconducting Josephson 
junction in a parallel magnetic field [3]. In co-planar double- 



ring BECs, such fluxons appear due to a competition between 
trap rotation and coherent tunneling. While the former de- 
mands a tangential velocity Q.R different for the two rings, the 
latter favors identical tangential velocities. Spatially nonuni- 
form or time-dependent potential differences between the two 
rings as caused by gradient or curved magnetic, gravitational, 
or optical fields generate forces on the fluxons. We also con- 
sider the conditions under which quantum tunneling of flux- 
ons may be observed. Further theoretical and experimental 
study of rotational fluxons in BECs promises to shed new light 
on the behavior of collective excitations in interacting quasi- 
onedimensional systems IU411 and their macroscopic quantum 
properties J3]. 

Below we start by presenting the theoretical description of a 
tunnel-coupled co-planar double-ring system, which is based 
on the Gross-Pitaevskii (mean-field) equation with a radial 
double-well potential. Its solution provides the ground-state 
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FIG. 1: (Color online) Phase diagram of co-planar double-ring 
BECs. We plot the difference of angular-momentum expectation val- 
ues for condensate atoms in the outer and inner rings as a function 
of rotation frequency Q, and tunnel coupling J. Finite integer values 
observed at higher frequencies are associated with the presence of 
rotational fluxons. Results shown are obtained for a typical double- 
ring geometry with d — 0.36 and g/ Jo = 100. Parameters and units 
are defined in the text. 
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FIG. 2: (Color online) Single fluxon in the ground state of a double- 
ring BEC, signified by the step-wise spatial variation of the relative 
phase a for condensate fractions in the two rings. Also shown are 
partial condensate densities | %i>i/o | 2 for the inner/outer ring. (9 is the 
azimuth. Parameters are Q = 5.8f2o, J = l-9Jo, and those used in 
Fig.E) 

phase diagram as a function of trap-rotation frequency and 
tunnel-coupling strength, as shown in Fig.Q] For finite tunnel 
coupling and a slow rotation, the phase difference between 
the two partial condensates in the individual rings vanishes. 
However, beyond a critical value of the rotation frequency, 
a quantized relative-phase winding between the two rings is 
accommodated, corresponding to a single rotational fluxon 
whose phase and density profiles are illustrated in Fig. [2] The 



ihd t ipi/o = 



Here S = (E - E { )/2 and = (E + E$/2 in terms of 
the single- well bound-state energies Ey . Rotation around 
the trap axis with frequency fi is imposed by any (initial) 
anisotropy in the trapping potential. Using the the normal- 
ization condition Y^ a =ioI l^a^dO = 1 the nonlinear cou- 
pling energies are gu — ng^\ where n = N/(2nR) is an 
average linear particle density and g[^ is the effective one- 
dimensional coupling strength Jish . For convenience, we in- 
troduce the effective trap radius R = \/2R Ri/ \J R 2 + Rf 
and d = (R 2 — R 2 )/(R 2 + Rf ), which is a measure of the ra- 
dial wells' separation, as parameters instead of i?j/ . We have 
solved Eq. (fl]i using an FFT-based pseudospectral method 
with imaginary-time propagation II 1611 to find the ground states 
of double-ring BECs. For simplicity, we assumed g\ — g = 
g. To compensate a trivial energy shift between states in the 
inner and outer well due to finite rotation, we have set 5 to 
5* = Mn 2 R 2 d/[2{l - d 2 )] for Figs. 1 and 2. 

In the absence of interactions (i.e., g — 0), stationary so- 
lutions of Eq. (Q~|i can be labelled by the quantum number 



phase structure can be detected experimentally by interferom- 
etry. E.g. after switching off the double-ring trap, both BECs 
will overlap in expansion and interfere destructively (resulting 
in a density node) at the azimuthal position of the fluxon. At 
higher critical values of rotation frequency, the number of ro- 
tational fluxons successively increases and a one-dimensional 
fluxon lattice is formed. After explaining our microscopic 
model and discussing numerical results, we present results 
from an effective hydrodynamic theory for BEC phase and 
density variables that captures the numerically observed be- 
havior. A classical equation of motion for rotational fluxons 
can be derived, showing that time-dependent and/or spatially 
nonuniform potentials accelerate fluxons. We find an expres- 
sion for the fluxon's inertial mass and discuss the possibility 
of quantum effects exhibited by these collective degrees of 
freedom. 

Theoretical description of co-planar double-ring BECs. We 
consider a situation where a BEC of atoms having mass M 
is confined by an external magnetic and/or optical trapping 
potential to two concentric rings with, in general, different 
radii R\ and R a for the inner and outer ring, respectively. We 
assume that transverse excitations in the individual rings are 
frozen out, allowing for a purely one-dimensional description. 
In addition, the two rings are coupled linearly by tunneling 
through a barrier with an associated tunnel energy J. Follow- 
ing Ref. iKj we consider the coupled Gross-Pitaevskii equa- 
tions for the inner and outer ring wave functions ipi(8, t) and 
ipo{9, t), respectively, which are given by 



■0i/ o - J i/Vi • (1) 

I 

hm of the angular-momentum component L z = —ihdg along 
the symmetry axis of the trap. The condensate wave func- 
tions in the inner/outer ring will be given by ip\j (9) oc e lm9 , 
and the phase difference between condensate amplitudes in 
the two rings will vanish at every point 0. However, a fi- 
nite g introduces a mixing of amplitudes with different m 
values in the condensate wave function, enabling the appear- 
ance of nontrivial structure in the relative phase. To illustrate 
this point quantitatively, we calculated the difference of ex- 
pectation values of L z per particle in the outer and inner-ring 
condensate fractions, i.e., (AL Z ) = (L Z ) Q — (L z )- U where 
(L z )n = (V^i/ol ^ IV'i/o) / (^i/ol^i/o)- In Fig. Q] (AL z )/his 
plotted as a function of tunnel coupling J [measured in units 
of Jo = h 2 1 (2MR 2 )] and rotation frequency ft [measured in 
units of f2o = H/(2MR )], for a particular double-ring ge- 
ometry. Regions with finite integer (AL z )/h are observed, 
which correspond to ground states with (one or more) rota- 
tional fluxons present. A representative example for such a 
fluxon's relative-phase and partial-condensate density profiles 
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is shown in Fig. [2] 

Basic features of the phase diagram shown in Fig. Q] can 
be understood by a variational consideration that assumes 
(i) strong nonlinear coupling g such that both rings are pop- 
ulated with equal density, and (ii) the condensate wave func- 
tion in each ring to be given by an L z eigenstate, ipy^^id) = 
e™ 1 ' 1 '" 6 1 \[^k . The values of m\i are determined by a competi- 
tion between tunneling, which tends to enforce equal phase 
for condensate fractions in both the inner and outer ring 
(to; = m = to*), and rotation. The latter favors the two 
condensate fractions to have, in general, different angular mo- 
menta determined by the rotation frequency and the ring radii 
(TO i/o = m i/o = lat[MRy ft/h+ 1/2]). It is straightforward 
to derive the energy functional of the system, 



£[to ,toi] = — [j^ + -j^J — ^-(m +mi)-JS m<> , mi . 

(2) 

The condition £ [to , TOi] = £ [to*, to*] defines a critical value 
J cr = MVL 2 R 2 d 2 /[2(l - d 2 )]. For J > J cr , the state having 
to; = to = to* would be expected to be the ground state, 
corresponding to the black region in Fig. Q] In the opposite 
case, the phase gradient for partial-condensate wave functions 
in the two rings will be different, essentially realizing a vor- 
tex (or several vortices) in the phase difference between the 
two rings. Such a situation is signified by the brighter colored 
regions in Fig.Q] The variational estimate of J cr yields a rea- 
sonably accurate description of the actual phase boundaries 
seen in the numerically obtained phase diagram. 

Effective analytical theory of fluxon phase profile and dy- 
namics. To obtain a more detailed understanding of fluxons 
in coupled annular BECs, we consider the dynamics of their 
collective phase and density variables. This approach applies 
equally well to co-planar and vertically separated double- 
ring traps. Writing the partial-condensate wave functions as 
ipilo = IV'i/ol exp{i(/>;/ }, we define symmetric and antisym- 
metric combinations of their modulus and phase and express 
the Lagrangian of the double-ring system in terms of these 
new quantities. It is possible to derive a closed equation of 
motion for the phase difference cf>- d — tp — <\>\ that is accu- 
rate to first order in the typically small quantity En/g, where 
Er = h 2 / (2M R 2 ) is the scale of energy quantization on the 



J 



ring. Its lengthy analytical expression is omitted. 

In the stationary limit and to leading (zeroth) order in 
Er I g, we find 



(1 - d 2 )E R d 2 (t). d - 2 J sin0 a = 
which has a single-soliton (i.e., fluxon) solution [1 

k(9 - O 
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Here sn(u|fc) is a Jacobi elliptic function [ 18] whose parame- 
ter k is determined from the transcendental relation 

TTK = kK{k) , (5) 

involving the complete elliptic integral of the first kind, and 
k = y / 2J/[(l — d 2 )E R }. Hence, fluxons emerge as station- 
ary phase configurations, as seen in our numerical calcula- 
tions. The dimensionless parameter n = R/(^/l — d 2 X.j) can 
be interpreted as the ratio of the quadratic mean radius of the 
trap R/s/l-d? = yJ{R 2 + R 2 )/2 and the physical length 
scale of the fluxon Xj = h/ (2\/ M J), which is set by the 
tunnel coupling. 

To obtain a dynamical equation for a slowly moving fluxon, 
we insert the ansatz <j)- d (9, r) = 4^(9, Oo(t)) into the equation 
of motion for the phase difference. Here 9o (t) is the instan- 
taneous position of the fluxon. Straightforward algebraic ma- 
nipulation yields a Newton-like equation of motion: 



Mfi9 R. — Ffi 



(6) 



The fluxon's dynamical mass is Mg — 2^J Ep / gl^M with 
the dimensionless moment of inertia given by 



2. = ii+d 2 )jj^ [*<ei 

_ 1 + d 2 K / 2-KK 

7r k \ k 



(7a) 
(7b) 



where E(u\k) is the incomplete elliptic integral of the second 
kind 1U8I1 . The general expression for the force (torque) on the 
fluxon is HI 




ddtf3 + d 



(1 - d 2 )E 



r 



2h 



(8) 



Equations dTbl and (0) define a universal relationship be- 
tween the fluxon's dimensionless moment of inertia Tn and 
the variable k. The limiting value of Tn for small trap size 
(k <C 1) is a constant [(1 + d 2 )/2], whereas a linear depen- 
dence [2(1 + d 2 )n/TT] is realized for large ring traps (k ^> 1). 



Inspection of Eq. ([H) reveals that fluxons subject to spa- 
tially nonuniform S and/or time-dependent S or (3 will expe- 
rience a force. This feature is confirmed by our numerical 
solution of Eq. (HJ, an example being shown in Fig. [3] In 
the case of spatially uniform 6(6, t) = So(t) and (3 = 0, the 
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FIG. 3: (Color online) Dynamics of a single fluxon under a time- 
dependent external field according to Eq. Q} (color coded density 
V ; i| 2 + V ; o| 2 )- An external potential gradient across the double-ring 
trap corresponding to S — 0.2/3 = 0.2 Jo cos(6>) is smoothly turned 
on at t = lOfi^ 1 and off again at t = SOfJ^ 1 . Other parameters are 
J — Jo and g = 100 Jo- 



Discussion and conclusions. We have discovered fluxon- 
like topological structure in the relative phase of condensate 
fractions in the ground state of BECs in rotating double-ring 
traps. These rotational fluxons are accelerated by spatially 
varying external potentials that couple asymmetrically to the 
two rings and/or time-dependent potentials. Macroscopic 
quantum tunneling of fluxons may become observable and 
would serve the long sought goal of preparing macroscopic 
quantum superposition states of BECs (see e.g. Ref. 12111 '). Fu- 
ture studies will focus on details of fluxons' quantum proper- 
ties and possible applications j^Q]. 

JB is supported by the Marsden Fund Council (contract 
MAU0706) from Government funding, administered by the 
Royal Society of New Zealand. 



force simplifies to y/2M/ gaS^, which is similar to the re- 
sult found previously 111 211 for phase-imprinted fluxons in a 
junction between two parallel linear BECs. The sign a = 



sgn 



^(27r)-$»(0) 



is the topological charge of the fluxon 
related to its orientation. Here we found the expression for the 
force felt by fluxons in the more general case with d 7^ 0. 

If the external fields are time-independent and the fluxon 
length Xj is smaller than the length scale of spatial varia- 
tions of 6, Eq. ((HJl can be integrated and written as Fy = 
—R~ 1 dgV, where V is a potential energy. For it > 1 and 
to leading order in d, we obtain 



V(6) 



7^g 




(9) 



for the potential and Ma = M ^32J/(n 2 g) for the dynamical 
fluxon mass. 

Macroscopic quantum tunneling. Describing the effects of 
quantum and thermal fluctuations on the fluxon dynamics can 
proceed in analogy to the established treatment of Joseph- 
son vortices in superconducting junctions y|. In particu- 
lar, the possibility of fluxon (macroscopic quantum) tunnel- 
ing can be included [20] by direct quantization of the clas- 
sical equation of motion Q6). A rough estimate for tunnel- 
ing of a fluxon through a potential barrier of height AV 
and length Al from the WKB method yields the probability 
P w exp(-2AV2A/ fl AV7ft). In order to have P > 1/e 
with Al w A j, we need \fJ~g > AV. Assuming a double-ring 
configuration as proposed in Ref. J3l with R ss 50 /im, it may 
be feasible to achieve g/ks ~ 2 /iK and J/ks ~ 0.05 //K 
and observe quantum tunneling through barriers AV/ks < 
0.3 /iK at sufficiently low temperatures. 

Quasiparticle excitations present at finite temperature will 
act as a damping mechanism for fluxon motion [3] and, at 
the same time, as a source of quantum decoherence (thus sup- 
pressing fluxon tunneling J20J]). 
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